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In this work, we evaluate the contributions of the hadronic loops to the amplitudes of J/i/; PV 
where P and V denote light pseudoscalar and vector mesons respectively. By fitting data of two 
well measured channels of J/t/j PV , we obtain the contribution from the pure OZI process to the 
amplitude which is expressed by a phenomenological quantity \Qs^ |, and a parameter a existing in 
the calculations of the contribution of hadronic loops. In terms of a and ISs^l, we calculate the 
branching ratios of other channels and get results which are reasonably consistent with data. Our 
results show that the contributions from the hadronic loops are of the same order of magnitude 
as that from the OZI processes and the interference between the two contributions are destructive. 
The picture can be applied to study other channels such as PP or VV of decays of J/^. 

PACS numbers: 13.20.Gd 



I. INTRODUCTION 

By the commonly accepted point of view, narrowness of the J/i/j resonance is interpreted by the OZI ruleP|. 
Namely, J /tp dissolves into three gluons which would then hadronize into measurable hadrons and these processes are 
OZI suppressed. In Fig.l, the threc-gluon process diagrams are shown and they are single-OZI-suppressed (SOZI) or 
double-OZI (DOZI) suppressed. Generally one can ignore the contributions from DOZI. 




(a) (b) 



FIG. 1: (a) and (b) correspond to the SOZI and DOZI processes for decays of J/t/j — > two light mesons. 

Due to complexity of the loop calculations and non-perturbative QCD effects which govern the hadronization, 
accurate estimate of the contributions from such OZI-suppressed processes to the decay rates is still missing. Therefore 
one cannot rule out other possible mechanisms which also make substantial contributions to the decay rates of J /ip 
and other family members of charmonia. 

Meanwhile, in the charm-tau energy regions, some phenomena cannot be understood in the regular theoretical 
framework. The pn puzzle is the most challenging one. Rosner recently suggested that this puzzle may be 
explained by the 2S — ID mixing. There exists an alternative possibility which might explain the surprisingly small 
branching ratio of -0(3886) pir, i.e. the final state interaction (FSI) 

Of course the first challenge to theorists is to correctly evaluate the widths of the exclusive J/tp non-lcptonic 
decays. Besides the OZI suppressed mechanism, there exists another possible process which may also contribute to 
the amplitude. That is the contribution from the hadronic loops. In fact, except the phase space for various channels, 
the mechanism where three gluons are exchanged, should be flavor-blind and universal for all the processes. It is 
another reason to motivate us to seek for a new mechanism because a simple analysis indicates that data of Jf-ip 
decays do not manifest the expected universality. 

This picture looks similar to that for the final state interaction (FSI) 0, Q, but essentially different. In the 
FSI picture, J/ip or -0(3770) decay into two real hadrons and the two real hadrons re-scatter into another pair of 
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hadrons which have the same isospin structure, but different identities, via strong interaction. In that case, the two 
hadrons, no matter in the intermediate stage or in the final states, are real and on their mass shells, therefore the two 
light hadrons cannot be charmed mesons such as D or D* due to constraints of momentum-energy conservation i.e. 
Mjj^ < 2M£)(or Md + Mo*)- One can derive the re-scattering amplitude by calculating the absorptive (imaginary) 
part of the so-called triangle diagrams 0, 0, Hi ■ 

There exists a different contribution. Recently, Suzuki indicates that the real part of the triangle where the two 
intermediate hadrons in the triangle can be off-shell virtual ones, does also contribute. Thus to estimate the corre- 
sponding contributions, one needs to calculate the dispersive (real) part of the triangle. The diagram corresponding 
to the hadronic loops is shown in Fig. 2 (there may be more similar diagrams which contribute to the processes (see 
Fig.3)). 




FIG. 2: The hadronic loop diagram where the two hadrons inside the triangle are time-like, but virtual ones and the exchanged 
hadron possessing appropriate quantum numbers are also off-shell. 

In fact, such a mechanism should also exist in all the decay modes of J/i]: and make sizable contributions. In 
this work, we are going to evaluate the contribution from the hadronic loops where the intermediate hadrons are 
off-mass-shell. 

The physical picture is following: J /tp first dissolves into two virtual hadrons, then by exchanging an appropriate 
hadron (i.e it possesses appropriate charge, flavor, spin and isospin) they turn into two on-shell real light hadrons 
which are to be seen by detector. The two intermediate hadrons are charmed hadrons, which contains flavor c or 
c, thus the transition J/V' H1H2 where H_\ and H2 are the two charmed hadrons, does not suffer from the OZI 
suppression. It is noted that the authors of discuss a possible mechanism for the relativistic heavy ion collision, 
where J /ip absorbs a pion existing in the hot environment to transit into two charmed mesons and the process may 
influence the observed J/V' suppression at RHIC. Even though such process is different from that under discussion of 
this work, the pictures have similarities. 

In the practical calculations, the coupling constant of the effective vertex cannot be obtained by fitting data because 
such a process cannot occur due to the limited phase space. Then certain reasonable symmetries are invoked to get 
the coupling constants. We can argue that even though the coupling constant itself is not accurately determined, the 
qualitative conclusion is not seriously infiuenced. 

In contrast with the derivation of the absorptive part of the triangle, the loop integration seems to be ultraviolet- 
divergent. In fact, following the standard procedure, the effective vertices are obtained from the chiral Lagrangian 
where all the concerned hadrons are on-shell, i.e. are real ones. To compensate the off-shell effects of the hadrons in 
the triangle, one may introduce form factors at each vertex. The form factors are inversely proportional to (g^ — A^) 
where A is a phenomcnological parameter Q which is expressed in terms of another phcnomenological parameter a 
(see the following text for details). Because of existence of the form factors the ultraviolet divergence disappears, 
namely A indeed plays a role similar to the cut-offs in the Pauli- Villas rcnormalization scheme |^ [lol | . 

Picture is as follows, because, so far, we cannot derive the contribution from the OZI suppressed processes (i.e. via 
the three-gluon intermediate stage) to the total amplitude based on the quantum field theory yet, we denote such 
contribution by a phcnomenological parameter l^/f^l which should be fixed by fitting data. 

Therefore, by the simple picture, there are only two parameters, i.e. a and |tj_g^|. Concretely, we suppose that 
both the direct three-gluon process and hadronic loops contribute to the amplitude altogether, by fitting data of two 
distinict channels of J PV where P and V stand for pseudoscalar and vector mesons respectively, we obtain 
the values of a and |t/5^| simultaneously. Here we choose to study processes J/ij] ^ PV, since there are more data 
available. Of course when we derive I^Jg^l from data, we have to carefully take care of the phase space of final states. 
Applying the obtained a and we evaluate the rates of other J/ip PV channels. 

Indeed, our numerical results show that a destructive interference between the three-gluon process(OZI) and the 
hadronic loop can result in rates which are satisfactorily consistent with data. We will give a detailed discussion in 
the last section. 

This work is organized as follows, after the introduction, in Sect. II, we formulate the hadronic loop contribution 
for J/ip PV. In Sect. Ill, we present our numerical results along with all the input parameters. Finally, Sect. IV 
is devoted to discussion and conclusion. 
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II. FORMULATION 

The effective Lagrangian for J /ip decaying into light pseudoscalar and vector mesons is written as 



(1) 



where G^^ , and G^^ respectively denote the effective couplings corresponding to single OZI diagram (Fig. 
(a)), double OZI diagram (Fig. ^ (b)) and hadronic loop diagram (Fig. In this work, we neglect the DOZI 
contribution. P and V respectively denote the nonet pseudoscalar and the nonet vector meson matrices 



K+ 



V = 



V K*- 



p+ K*+ 



(2) 



(3) 



In the above expression, Xrj^r]' and y,;,,/ describe the rj — -q' s mixing, and 



cos 9 — \f2 sin Q 



sin 6* + \/2 COS0 



(4) 



where B ~ — 19.1°[T3| is the mixing angle of r\ and 77'. 

Taking J/?/) — > (P'k^ as an example, we illustrate the calculations of the contributions from hadronic loops. 

The Feynman diagrams describing the contributions of the hadronic loops to J/i/' ^ p°7r° are depicted in Fig. O 
In the hadronic loops of diagrams (a)-(f) only and Z)(*)o are involved, but of course they can be simply replaced 

by to compose new diagrams. 




{a) 




(6) 




J/V' 




{d) 



(e) 



(/) 



FIG. 3: The Feynman diagrams which depict the decays oi J /ip —>■ p'^n'^ through D'*' hadronic loops. 
In the former works the effective Lagrangians, which are related to our calculation, arc constructed based 
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on the chiral symmetry and heavy quark symmetry as 

+tg^,^,,V*^^d^.V*i{Vn - + 4z/,.,.,P%ya^V'' - d-'Wy^vA, (5) 



where V and V* are respectively pseudoscalar and vector heavy mesons, i.e. V^*^ ^{{D^)'^*\ {D~Y*\ {DjY*^). The 
actual values of the coupling constants will be given in the following section. 

With these Lagrangians, we write out the decay amplitudes of J/t/j p°7r°. For Fig. El (a), 

i i(—q°"^) i 9, 
Pi -iTT-oPi- '^D' 9 - ™ri 

For saving the space, we collect the concrete expressions of the amplitudes corresponding to Fig. (b)-(f) in the 
Appendix. 

In the amplitudes, Apo^^o = 1/2. J- {mi, q^)s represent the form factors which compensate the off-shell effects of the 
mesons at the vertices and are written as [l5j 

9n — mf 

where A is a phenomenological parameter. It is obvious that as ^ it becomes a number and if A ^ m^, it turns 
to be unity. Whereas, as — *■ 00, the form factor approaches to zero. In that situation, the distance becomes very 
small, so that the inner structure would manifest itself and the whole picture of hadron interaction is no longer valid. 
At that region the form factor turns to zero and plays a role to cut off the end effect. The concrete expression of A 
is represented as 

A(to,;) ^ rrii + aAgcD, (8) 

where rrii denotes the mass of exchanged meson. Aqcd = 220 MeV. a is a phenomenological parameter. In fact, in 
literature, some other forms for the form factors are suggested fl^. all of them depend on phenomenological parameters 
which are obtained by fitting data, therefore in the calculations all the forms are somehow equivalent even though 
some of them may provide a better convergent behavior for the triangle loop integration. That is very similar to the 
case of the pauli- Villas renormalization scheme 0, ■ 

Taking Fig. |31 (a) as an example, we carry out the Feynman parameter integration for © and get 

' ''^^ ' 2 ^^^^A(™„,m„.,A(™J) 



where 



(A2 - ml)y 



87r2A(mo,m^. ,TOo) 



(47r)2 V A(mj3, m^^, , m^^) 
(-2£^.a;3e^/^e>^pf), (9) 



A(a,5, c) = 7713(1 — .T — y)^ — 2(p3 • p4)(l — X — ?;).T + 77142;^ — (mg — a^)(l — X — y) 
-(to| - b'^)x + yc^, 

nT-sin^i), PsiPA) are the masses and four- momenta of p^{tt'^) respectively. Because of the form factors, the ultraviolet 
divergence does not exist as expected. 

We use the same method to deal with the amplitudes corresponding to Fig. O (b)-(f) and also put them in the 
appendix. 
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To obtain the amplitudes corresponding to hadronic loops which contain D^*')^, one only needs to replace the 
parameters corresponding to and Z)(*)° by that to and D^*^~ in the above six expressions l|^- H27|) . Thus 

contribution from the hadronic loops to the total amplitude can be eventually expressed as 



Mi^" = {Ma + Mi, + M, + Md + Me + Mf) 

+ {Ma +Mb+Mc + Md + Me+M f)\r,^^^,^o-.m^M± 

- Gi'^\^)e,.c.pe':j^^e;p^pi (10) 

where Qi^^° [a) denotes the effective coupling for J/ip ^ p°7r° induced by the hadronic loops where Q'^^" {a) is a 
function of variable a. 

Similar to J/tp — > p°7r°, we can also obtain contributions of hadronic loops to J/ip — > K*~^K^ + c.c, K*'^K'^ + 
C.C., Lurj^ \ (pT]'- \ It is noticed that the hadronic loops for processes J/ip — > K*'^ +c.c., K*°K^ + c.c., (jyrj^ ^ include 
ui*^^^^ mesons and in the calculations, the coefficient ApO^o is replaced by Ak'K ~ 1 (or ^^^o = A^^^o = 

We have obtained seven quantities which contain both contributions of the OZI process and the hadronic loop 

G""-" - gri3'^"''° + g'H^\a), (11) 

g-" = e|'^/?'"'' + g^'" + g7(a), (13) 

= (14) 
S"^" = 5^13"^'^ + gp + gp{a), (15) 

g't'v' ^ g^^'is^'"' + gf,''' + gp'ia), (i6) 

where the contributions from the DOZI processes are ignored, i.e. one can set g'^^ ' = gf^^ ' = 0. Geometry factors 
(3^^ satisfy the following relation: 

/3f : Pf^""' : ■■ ■■ ■■ = 1:1::^.:^.': : y,' ■ (17) 

The SU(3) geometry factors /3^^ is obtained by taking the trace Tr{PV) as J/i/' is an SU(3) singlet, and in fact, it 
is just the SU(3) C-G coefficient for the concerned channel. 

Thus using eqs. (|ll(l and (|12|l . the parameter a is determined. With it, we are able to calculate the values of g^ ■ 



III. NUMERICAL RESULTS 

Because of the flavor SU(3) symmetry, it is reasonable to set g,/^j^ ^ ~ 9 1/^,00- Based on the naive vector 
dominance model and using the data of the branching ratio of J/tp — > e^e~, the authors of Ref. [l^ determined 
(7^^^^^/(47r) = 5. As a consequence of the spin symmetry in the HQET, ,9^^^^ ^, and g j/^j^,^^ satisfy the relations: 
9j/^n,nt = 9j,^nnJ^n, and = 9.,,^n,n, El- O^^^cr Coupling constants related to J/i/' PV include [3: 



~~ 9d*Dt,1 9d*DK " \ ZZ 9d*D7,^ 9v ~ f 1 

where /^r = 132 McV, (3 and A are parameters in the effective chiral Lagrangian that describes the interaction of 
heavy mesons with low-momentum light vector mesons 0|. Following ref. (23, we take g = 0.59, /? = 0.9 and A = 0.56 
GeV^ll^. Meanwhile, the coupling constants .) (/n(')£,(*)0(^.)) are related to g dm Di') pif Di') dm p) and 

one can expect the relations in the limit of SU(3) symmetry [ig 

= —^9DMDMp. fDi''Di'>4> = ZT^foMoMp, 



TDMn^'U<* ^ 1/ ~gDMDMpi fnMnMf^,—\ —foMDMp- 
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Using cqs. and H12() . we have determined a and as 

a ==0.13, l^^l"^! = 4.51. 

Namely, in our scheme, the hadronic loop contribution to J/ij] ^ PV is theoretically calculated whereas the OZI 
contribution is obtained by fitting the data of the two channels. The obtained \Gs^ \ = 4.51 is supposed to be universal 
for all channels of J/V' — *■ PV based on SU(3) symmetry. With eqs. H13|l - H16() . we predict the values of Q^^ for all 
channels oi J/ip ^ PV. The obtained results are listed in Table. ^ 



Decay mode 




K*+K- +C.C. 


4>V 


4>v' 


LUT] 


cur]' 


BRx 10"^ (Experiment) [21] 


4.2 ±0.5 


5.0 ±0.4 


0.65 ± 0.07 


0.33 ±0.04 


1.58 ±0.16 


0.167 ±0.025 


5^^ (10-=^ GeV-i) 


2.08 ±0.25 


1.65 ±0.26 


0.89 ± 0.096 


0.71 ± 0.086 


1.27 ±0.13 


0.46 ± 0.069 


g^^" ( 10"^ GeV" ' ) (Theory) 


6.59 


6.16 


3.55 


5.13 


5.46 


4.07 


g^^(10"^GeV-^)(Theory) 


2.08(fitting) 


1.65(fitting) 


0.92 


0.62 


0.95 


0.40 


BRx 10"^ (Theory) 


4.2(fitting) 


5.0(fitting) 


0.70 


0.25 


0.86 


0.14 



TABLE I; The first two modes are well measured, we obtain a and [Gs^' \ = 4.51 by fitting them. 



IV. CONCLUSION AND DISCUSSION 

From above discussions and numerical results, one has the following observations. 

(1) The contribution from the hadronic loops is sizable and has the same order of magnitude as that from the SOZI 
process. 

(2) To get reasonable results which are consistent with data, the interference between the SOZI contribution and 
that from hadronic loops is destructive. 

(3) We derive the parameter a and the universal l^/g^l by fitting data of the first two distinct modes of J/tp and 
present in the first two columns of table 1. In our picture, there are only two free parameters. The values of 
which we predict in terms of a and Gjp^ are reasonably consistent with the experimental data. 

In this work, we rely on two key assumptions. The first is that the value l^/f^l is universal for all channels of 
J/^p — > PV due to the fiavor-blindness of gluon coupling to quarks. This treatment is similar to that made for 
study process tt + J/ip Z?^*-* + -D'*-* in a hot RHIC environment. 

In this work, we re-parameterize A in terms of a following the ansatz suggested by the authors of Ref. [l3|. It seems 
that one may gain some information about the parameter a from the FSI calculations 0,1^, however, it is not really 
the case. Even though in both calculations, one deals with the triangle diagrams which are composed of only hadrons, 
the situations are different. For FSI, one calculates the absorptive part where only the exchanged meson is off-shell, 
whereas for the calculation of the dispersive part, all the three mesons in the triangle are off-shell. Form factors are 
phenomenologically introduced to compensate the off-shell effects at each vertex, thus the form factors should not 
be the same for calculating the dispersive and absorptive parts, i.e. even though they have the same expressions, 
the involved free parameters may be different. It is noted that in the FSI case, the value of a in the form factor is 
expected to be of order unity 0|, but in the case of hadronic loop the value of a is obviously smaller. In our strategy, 
we do not priori set its value, but let data determine it. 

As our conclusion, it is obvious that the contribution from the hadronic loop to J/tJj two light mesons should 
exist. The question is how large it is compared to that of the SOZI processes. Namely, if only the SOZI process is 
responsible for the decays of J/tp, it is hard to explain the data and as the contribution from the hadronic loops being 
added, one can see that the data are well accommodated within a tolerable range. In our scenario, there are only two 
parameters in the whole scenario and they are determined by fitting data of two channels. With these two parameters, 
we evaluate the branching ratios of the other channels oi J/tp ^ PV and the theoretical predictions are consistent 
with data. It indicates that the ansatz is reasonable and the whole picture makes sense. Indeed, there is only one 
channel, i.e. J/tp tor], for which the theoretical prediction on the branching ratio is about 0.6 of the experimental 
data. Even though there exists such a discrepancy, for a rough estimate the consistency is not bad at all and we 
would like to urge our experimental colleagues to redo the measurement to get more accurate results. On the other 
side, so far the SOZI processes have not well calculated in the framework of quantum field theory yet, we also need to 
complete such a derivation and detect the value which we obtained by fitting data and assuming a universality for all 
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modes. To be more accurate, however, further theoretical and experimental efforts may be needed. This mechanism 
should exist not only for J /ip exclusive decays, but also for other members of charmonia. However, for higher excited 
states of chamonia, the channels with final states composed of charmed hadrons are open, and the situation becomes 
more complicated. This picture also applies to the T decays and can be tested by the future experimental data. 
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Appendix 



One obtains the decay amplitudes of J/ip — > p°7r° corresponding to Fig. |2l(b)-(f). 
For Fig. El (b). 



Mb = Aa 



P«^n j j2:^^'^3J/i-DD-£t.^ai3e'}^^{ip2){iPi)][-2ifD'-Dpe'''^ ^ {ip^)el{-ip^ - iq )] 

^[igD'D*TT£7TCSLj{ipi){~iP2)]^ T 2 ^ 2 9^ ^ 2 ^^{mD',q^)- (18) 

Pi ~rnjjP2- ^D- 1 - mjj. 



ForFig.El(c), 



j^[i9j/i,DD{pi-P2) ■ ej/^][-2ifD'Dpe''^^^{ip^)el{~ip^ - iq^)] 

^[gO'D^itp^)] o ^ 2 9^9 'i_l72^ T\mD',q^). (19) 



For Fig. 01(d), 



Pi - mj^ P2 ~ "^b 9 ~ "^b* 



J^^iaJli'D-D' [iPi^^j/^, - V2^j/i, - i{P2 - Pi) ■ 
^[2ifD- Dp£iTCSi^{ipl)e},{-ipi - iq^)][gD' DT:{ipt)] 

X— 2 2 2 [mD,q . 20) 



ra. 



For Fig. El (e). 



For Fig. El (f), 



-J^^{9.JI^D'dA^p\('^jI^ - V2<^j/4, - «(P2 - Pi) ■ ej/.^S^"]} 
x{-.gD-_D>(-jpi - iq) ■ tpg""^ - ^fD'D'p{ipt4 ~ ^Pf^p)} 

AWD'D'.e^CS^{zpi){-ipi)] ^"^7 2^~^^2 P\^U'.q^). (21) 

Pi - mf). pi - mf). q^ - mf^. 



j^[igj/4'DD'£^L,yapej/^{iP2){iPi)] 
x{~9D*D'p{-ipi - iq) ■ epg^'' - ifo'D-piipUl ~ ip^e^)} 

X [9D'D. {^pt)] i^D-.q^). (22) 
Pi - m^y, pi -m^jjq^- mjj. 
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By further calculations, we obtain 



Jo J a 
1 1 



(A2-m|5.)2/ 



167r^A2(TOj-, , m^, , K{'m^, )) 



167r2 A(to^ , m^, , m^, ) 167r2 A(m^ , m^. , A(m^. )) 
'^x^yipz ■ Pif + '2xy{p3 ■ Pi)^ + 2xy'^m\m'l + 2x^1/77137714 - 2x2/777.37714 
2 . /A(77i^,77i„.,A(m^.))N^ (A2-m2,,)y 



^ 1 (' 



{-6y -2x + l){p3 ■ Pi) + 2ml - '^mlx - 2mly + 2y{p3 ■ p^) 



(23) 



{A^-ml,)y 



1 i-l-x 

dx / dy 
Jo 



2 



A(771„,777^,A(7?7^.)) 
A(r?7^, 771^, 777^.) 



87r2A(77lo, 771^, 777^.) 



(24) 



1 



1 l-l-X 

dx / dy 
Jo 



(A2 - ml))y 



1 



167r2 A(771^. , 777^. , 777^) 167r2A(77l^. , 777^. , A(7?7^ )) 
\2 S 



+4:xy{p3 ■ Pi) 



12 



(47r)2 



167r2 A2(777^. , 777^. , A(777o )) 

4x^7774 + 4a;^(p3 • pi) 

(A2 - 7772,)y 



A(777^. , 777^., 777^) / StT^ A(777^. , 777^. , 777^ ) 



(25) 



M{e) = ^p07r°9j/^o'n'9o'n'. I dx j dy 
1 



(A2- 7772,.)y 



1 



f . . 



IGn^Alm^, , 777 j-,, , 777 j-,, ) 167r2A(777j-,, , 777^, , A(777^, )) 
+4qj^m\ + 47^73 — 8x'7773 — 87/7(73 + Sxy7n\ + 4x^(^3 • p^) + 4(p3 • p^) 



-8x(j>3 ■ Pi) - 4y{p3 ■ Pi) + Axyips ■ pi) 

(A2 - 777|,.)y 



87r2A(777^. ,777^. ,777^. 

1 



(47r)2 V A(777^. ,777^. ,777^. ) 



1 j-l-X 

no*. I dx / dy 
Jo 



(A2- 777,2,, )y 



1 



167r2 A(7n^. , 777^ , 777^. ) 167r2A(777^. , , A(77i^. )) 



167r2A2(7nj^, , 777 j-, , A(7r7j^, )) 
- Ay{p3 ■ Pi) 



(26) 



(27) 
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